Comment on: 


"High-Pressure Elasticity of a-Quartz: Instability and Ferroelastic 

Transition" 
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Gregoryanz et al. ||l| have analyzed the stabihty of a- 
quartz under pressure in terms of the Born stability cri- 
terion (BSC) 1^. The BSC requires a positive definite 
elastic matrix constant C^ for stable elastic solids. It is 
strictly valid only at zero external pressure. Important 
differences were discussed with theoretical predictions by 
Binggeli and Chelikowsky 0: The pressure regime in 
which a-quartz could be expected to be mechanically sta- 
ble if it was possible to suppress the quartz I - quartz II 
transition at p = 22 GPa was found to be significantly 
larger than predicted theoretically. This conclusion is 
incorrect, because Gregoryanz et al. use an improper 
generalization of BSC to non-zero pressures. The correct 
criterion is to require that the matrix of Birch or stiff- 
ness coefficients [0 Bij is positive definite rather than 
Cij H. Binggeli and Chelikowsky also state their BSC 
in terms of Cij, but as it seems to us, they accidentally 
calculated Bij, e.g., they might have inverted the matrix 
V{6ei6ej)/kBT with Sci being the thermal fluctuation of 
the strain tensor around the equilibrium state at a given 
externally applied stress. This inversion gives the Birch 
coefficients rather than Cij; the equality Bij = Cij only 
holds at zero external stress. As we will show now by un- 
published molecular dynamics simulations of the present 
authors, the predictions of Binggeli and Chelikowsky are 
correct. We use the potential energy surface from Beest 
et al. [^, which is known to be particularly well suited 
to describe (elastic) properties of a-quartz. We focus on 
the BSC 


Bz{Cij) = (Cii — Ci2)C44 — 2C;^4 > 0, 


(1) 


which is shown in Fig. 1. 

In the inset of Fig. 1 it can be seen that the present sim- 
ulations reproduce the experimentally measured pressure 
dependence of B^[Cij) (the central Fig. 3c of Ref. [1]) 
fairly well. We also want to note that we obtain the 
quartz I quartz II transition at the same pressure p ± 
0.5 GPa as in experiment. These agreements justify the 
use of the BKS potential for a-quartz under pressure. If 
Bafj replace C^p in Eq. y, the behavior is similar to the 
one reported by Binggeli and Chelikowsky |^, e.g., a 
nearly vanishing slope of B^(Bij) at small pressures is 
found, and instability of a quartz can be expected to oc- 
cur at p « 25 GPa. This is the correct generalization of 
BSC, which is equivalent to require that all eigenvalues 
of Bij are positive. The minimum eigenvalue of Bij as 
obtained in our simulations is shown in the main part of 
Fig. 1. 
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FIG. 1. Minimum eigenvalue of the matrix of the elastic 
constant (open diamonds) and of the Birch matrix (filled cir- 
cles) as a function of pressure. Inset: Born stability criterion 
applied to elastic constants (diamonds) and Birch coefficients 
(closed circles). Error bars are about symbol size. Crosses 
represent experimental data from Ref. [1]. 


We further want to point out that the observation of 
dCii/dp < 0, which is discussed exphcitly in Ref. ||], has 
no direct implication on the mechanical stability. The 
only relevant influence of C44 on mechanical stability in 
the case of an isotropic pressure p enters via the properly 
generalized Born stabihty criterion B3{Bij), see Eq. (||), 
with B44 ~ C44 — p. 
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